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Abstract 

Near-critical binary mixtures containing ionic solutes near a charged wall preferentially adsorbing 
one component of the solvent are studied. Within the Landau-Ginzburg approach extended to 
include electrostatic interactions and the chemical preference of ions for one component of the 
solvent, we obtain a simple form for the leading-order correction to the Debye-Hiickel theory result 
for the charge density profile. Our result shows that critical adsorption influences significantly 
distribution of ions near the wall. This effect may have important implications for the screening 
of electrostatic interactions between charged surfaces immersed in binary near-critical solvents. 
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I. INTRODUCTION 



A phenomenon of critical adsorption occurs when a fluid is brought to its bulk critical 
point in the presence of an attractive substrate or wall, for example, along the critical 
isochore. The wall causes a perturbation of the relevant order parameter profile to extend 
over a distance ~ the bulk correlation length, from the surface l|, |2[. Close to criticality, 
where ~| (T — T c ) /T c \~ u [y is the critical exponent), the influence of the wall extends 
to macroscopic distances. As a result the amount adsorbed (adsorption T) diverges as r = 
(T — T c )/T c — >■ 0. Here we focus on the equivalent phenomenon that occurs for binary liquid 
mixtures near their consolute points. In these systems, generically there is a preferential 
adsorption of one component of the mixture on the wall surfaces, resulting in the divergence 
of the relative adsorption near the critical point. 



and experimentally. Of interest 
that T should have a power law 



Critical adsorption was much studied both theoretical . 
has been the scaling prediction by Fisher and de Gennes I 
dependence on r, with a universal exponent which does not depend on the details of the 
specific system. Also, this phenomenon is of significant practical importance, e.g., for the use 



of supercritical fluids in micro- and nanofluidics 

0. 



4| and as solvents in colloidal suspensions 



In colloidal suspensions one usually encounters the presence of ions. This may be due 
to the dissociation of a salt-free solvent, e.g., in binary solvents containing water, or due 
to dissolved salt. Ions influence the electrostatic interactions between the charged colloids 
leading to screening effects If there is a chemical preference of ions for one component of 
the solvent, then larger amount of ions may be dissolved in regions where the concentration of 
the preferred component of the solvent is larger. For this reason critical adsorption in binary 
solvents may affect the distribution of ions near charged surfaces. This effect in turn may 
influence the screening of electrostatic interactions between charged colloidal particles. Such 
a scenario is relevant for an experimental system used recently to investigate the effective 
forces acting on spherical particles close to a substrate immersed in a near-critical water- 
lutidine mixture 3, |8|. The aim of that study was to measure directly a critical Casimir force 
between a single spherical colloid and a flat surface. Polystyrene particles with high surface 
charge density were used as colloids, and silica glass treated chemically to achieve the desired 
adsorption properties, which also carried some surface charge, was used as the substrate. 
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Ions in solution were present due to dissociation in the salt-free water-lutidine mixture. 
In order to extract the critical Casimir force from the measured potentials of interactions, 
the distances z from the substrate for which the electrostatic contribution was estimated 
to be negligible have been considered. The estimation of the electrostatic interactions was 
based on measurements of the electrostatic potential far from the critical point. However, 
different solubilities of ions in water and lutidine and the enhancement of one of components 
of the near-critical mixture close to the surfaces might result in a change of the screening 
of the electrostatic interaction compared to the case of a noncritical homogeneous medium. 
Unfortunately, possible effects of the critical adsorption on the electrostatic contribution to 
the interaction potential could not be estimated due to the lack of any theoretical predictions 
for such effects. 

Motivated by these issues in the present paper we study critical adsorption in a simple 
generic system, i.e., for a binary mixture solvent containing ionic solutes in the presence of 
a single charged planar wall that preferentially adsorbs one component of the solvent. 

Recently, the properties of bulk near-critical binary mixtures with ions have been studied 
theoretically using a modified Poisson-Boltzmann theory by Onuki and Kitamura 9J. The 
Landau- Ginzburg functional was used to describe the binary solvent. The presence of ions 
was accounted for by including electrostatic contribution and the terms due to the entropy of 
ions. Preferential solvation effects were modeled by a linear coupling between the density of 
ions and the order parameter of near-critical binary mixture. Using this phenomenological 
functional in the mean-field approximation, the authors found that the presence of ions shifts 
the critical point of demixing. This result is consistent with the well known experimental 
observations that adding salt to the binary mixture solvent changes the position of its phase 



separation curve, such that the upper and lower critical points are shifted upwards 
and downwards 12| . respectively. 
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If the ion-containing mixtures are exposed to external electric fields, e.g., due to charged 



walls or chargec 
is also shifted 



macromolecules immersed in the solution, the phase separation temperature 



131 ] . In addition to the preferential solvation, a dielectric inhomogeneity 



is important for this effect; the high permittivity solvent component is attracted to the 
charged surface thus enhancing the phase separation. Recently, the effects of these, so 
called, dielectrophoretic forces and of the preferential solvation have been studied for ion- 
containing mixtures confined between two-charged surfaces in a single phase region away 
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from the coexistence [ijj]. The approach used for a bulk system in Ref. {9] have been 
modified to account for a spatial variation of the volume fractions and, consequently, of 
the dielectric permittivity of the binary mixture. In addition, a surface term describing the 
interactions between the charged solutes and confining charged surfaces has been included. 
Using a mean-field approximation the density profiles and osmotic pressure between charged 
interfaces have been calculated. 

In order to investigate the influence of critical adsorption of ion- containing binary mix- 
tures on the distribution of ions near a single charged wall, we propose here an approach 
that starts from the microscopic theory. Within such an approach the Landau functional is 
derived rather than postulated, and it should describe correctly all the collective phenomena 
in the system that are consistent with the assumed interaction potentials. Hence, we can 
control all the assumptions on the fundamental level of interactions. In Sec. Ill Al we intro- 
duce the lattice-gas model of a four-component mixture describing two species of the solvent 
and the positive and negative ions. The generic for this model, short-ranged interactions 
between all the species and between the species and the surface are assumed. From the 
lattice model we derive the continuum Landau-Ginzburg functional which is then supple- 
mented by electrostatic bulk and surface contributions. Contrary to previous approaches 
all parameters of our functional are expressed in terms of microscopic interactions. Our 
Landau functional is similar to the functionals in Refs. js), [l^, but more terms are present, 
and all the approximations are based on assumptions concerning types of interparticle in- 
teractions. In Sec. IIIBl we derive Euler-Lagrange (EL) equations that allow to calculate the 
density and charge profiles on a mean-field level. Next, in Sec. II III we solve the linearized 
EL equations (Sec. IIII Al) and calculate the leading-order correction to them (Sec. IIII Bp 
under the assumption of weak interactions with the surface, small amount of ions and weak 
surface charge. A discussion of our results is included. The summary and outlook are given 
in Sec. ED 
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II. DERIVATION OF THE LANDAU THEORY 



A. Construction of the lattice model 

Let us consider a two-component mixture approaching the upper critical point of the 
demixing transition from the one-phase side and investigate the effect of the presence of a 
small amount of positively and negatively charged ionic solutes. We are interested in such a 
system in contact with a charged wall, such that the total charge of the system neutralizes 
the charge at the wall. Let us first consider this system with the ions replaced by the 
corresponding neutral molecules. The electrostatic contribution to the free energy will be 
included in the next step. 

The four-component mixture close to the demixing transition can be conveniently studied 
in the framework of the lattice-gas model. In the simplest lattice-gas model of a one- 
component system the space is divided into cubic cells of a volume comparable to the volume 
of a molecule. The cells can be either occupied or unoccupied by a center of a single molecule, 
i.e., multiple cell occupancy is forbidden. For a mixture of molecules of different sizes we 
consider a semi-microscopic description of cells occupied either by the larger molecule or by 
a cluster of smaller molecules, and we denote the volume of the lattice cell of a simple cubic 
lattice (SC) by Vq. Here we limit ourselves to a mixture of molecules of similar sizes. The 
cell-occupancy operators are Oj(x) = 1 if the cell x is occupied by the i-th component and 
6j(x) = otherwise. The two components of the solvent are denoted by i = 1,2, and two 
additional species that can be ionized to carry positive or negative charges are indicated by 
i — 3,4, respectively. We restrict ourselves to thermodynamic conditions corresponding to 
the stability of a liquid, where total density fluctuations can be neglected. The single-state 
occupancy and close-packing lead to the constraint 

X>(x) = l. (1) 

i=i 

We consider an open system with the probability of the semi- microscopic state {6j(x)} given 
by 

p[{Si(x)}] = exp(-^[{6 8 (x)}]) n x ^(Eti o_M - 1) } (2) 

where (3 = (/cbT) -1 , ks is the Boltzmann constant and T the temperature. H is the par- 
tition function. We include the activity factors in the Boltzmann factor, by assuming the 



generalized Hamiltonian 

H SR [{d^)}] = £ SR [{o^)}] - /*i£>(x) - /i 2 ^6 2 (x) - ^5> 3 (x) + o 4 (x)), (3) 

XXX 

where ^ is the chemical potential of the species i. We shall require charge neutrality for 
ionized species 3 and 4 in the bulk, and assume /J3 = For the short-range interaction 
energy we make the standard nearest-neighbor (NN) approximation for all the species 

£ SR [{d^)}] = -^^(x^x + e^ (4) 

x n=l 

where e n is the unit lattice vector in the n-th direction and summation convention is assumed 
for i,j = 1,...,4. The coupling constants are symmetric, Jjj = Jji, and represent the sum 
of all short-range interactions. The Hamiltonian of the lattice model with the electrostatic 
interactions disregarded takes the form 



h sr [s^]=-y: 

<x,x') 



s(x) J ss <S(x') + 2 (x) J pp 2 (x') + 0(x) J^0(x') (5) 



+s(x) J sp 2 (x') + 0(x) J 0s s(x') + 0(x) J 0p 2 (x') 



- ^ I A^s(x) + A/i c 2 (x) I - I M(x) + h p <f) 2 (x) J + const, 
xev V / xeay V / 

where * s the summation over all pairs of nearest-neighbors belonging to the system. 

The system volume is denoted by V, dV denotes the boundary layer, and the new operators 
are 

s = oi - o 2 = — 1, 0, 1 (6) 
<t> = o z - 6 4 = -1,0, 1 

4> 2 = <3 3 + o 4 = 0, 1. 

From close packing ([!]) it follows that 61 + 62 = s 2 = 1 — <p 2 . The new coupling constants 
in (0) are linear combinations of the coupling constants Jy, and A/x and A/j c are linear 
combinations of 

We further assume that the species 3 and 4 are similar, i.e., interact in a similar way with 
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the solvent and with each other, such that we can assume 



<^33 + J44 ~ 


2 J34 




a J 14 


J23 ~ 





(7) 



With this assumption the Hamiltonian takes the simpler form 



H s,i \S,i 2 ] 



EE 

xGV n=l 
2, 



s(x)J ss s(x + e n l 



fx) J, 



/>,<>■ 



,x 



+ 



(8) 



s(x)J sp (f) (x + e„) + <p (x) J sp s(x + e 




Note that the SR contribution to the energy depends only on 2 when the assumptions ([7]) 
are valid; otherwise the full expression ([5]) has to be considered. We assume J ss ^> J pp , 
so that the system can phase separate into two phases, one of them rich in the first- and 
the other one rich in the second component of the solvent, rather than into a phase rich in 
the 3-rd component and a phase rich in the 4-th component of the mixture. The coupling 
J sp > signals that ions preferentially dissolve in the first component of the solvent. 

When the species 3 and 4 are charged, then Coulomb interactions between them are 
present in addition to the short range interactions. In the truly microscopic model one should 
take into account polarizability of the solvent, and the problem becomes very difficult. In 
the semi-microscopic description the effect of polarizability of the solvent is included only 
through the dielectric permittivity. The dielectric permittivity depends on the composition 
and, in the case of spatial inhomogeneities of the composition, the dielectric permittivity is 
a functional of the solvent densities and depends on x. The electrostatic energy- density e e 
in this case has the form {{], flo1 | 



e e [s(x),p c (x),0(x),^(x)] 



8tt 



(V^) 2 + e0(x)^(x)] 



(9) 



where 



s x 



(*(x)>, 



(x)> 



(10) 
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are the mean values of the microscopic operators, V is the gradient and i/j is the electrostatic 
potential that in equilibrium corresponds to the minimum of the energy (jH]), i.e., to gener- 



alization of the Poisson equation, considered for example in Ref. Close to the demixing 
point of the binary mixture the densities vary on the length scale of the bulk correlation 

1/3 

length £b ^> v , and their deviations from the average value are small. It is then reasonable 
to assume that e(x) is a linear function of the average densities of the two components of 
the solvent 

e(x) = eipi(x) + e 2 p 2 (x) (11) 

where Pl (x) = (ox(x)) = (1/2)(1 - p c (x) + s(x)) and p 2 (x) = (6 2 (x)> = (1/2)(1 - p c (x) - 
s(x)). ei and e 2 are the dielectric constants of the pure species 1 and 2 respectively. When 
the electrostatic energy has the above form, the average densities have to be determined 
selfconsistently. 

If one assumes, however, that e(x) can be approximated by its average value 

e(x)«e, (12) 

then the electrostatic energy (Q has the same form as in vacuum (with the modified per- 
mittivity) and can be written as a sum of interaction energies for all pairs of point charges, 

£ C [{o t (x)}] = \J2J2 £ ^(x)K(|x - x'|)eA(x') (13) 

x x'^x i J=3,4 

where e% and are the charges, and V c is the electrostatic interaction potential., that has the 
form different for different lattices and in the continuum space. The forms of V c for simple 
cubic, body centered, or face centered lattices in the Fourier representation can be found, 



for example, in Ref. [161. Il7l ]. If we assume for V c the continuum-space form and restrict the 
positions of ions to the lattice sites, then the electrostatic energy takes the familiar form of 
the sum of Coulomb interaction potential for all the pairs of ions 

£ c [W(x)}l = f££|rr^ d4) 

x x't^x 

In the above expression we limit ourselves to monovalent ions e^, = —e^ = e, with e denoting 
the elementary charge. With the above form of the electrostatic energy we obtain a well- 
defined semi-microscopic model with the probability distribution 

^ = !^(zMMim^(f!±i!^), (15) 



s 



where "H[{s,0}] = H SR [{s, <f> 2 }] + The grand potential in the above statistical- 

mechanical model is given by 

tt = -k B T\nZ. (16) 

The advantage of such kind of modeling with the assumption ffT2l) is the possibility for 
investigating the effect of fluctuations of the composition (for example by means of computer 
simulations), and the disadvantage is neglecting the coupling between e (and hence the 
electrostatic energy of states {oj}) and the composition fluctuations. Which effect, the 
fluctuations of the composition and their coupling to the density of ions or the spatial 
variation of e, plays the dominant role in determining the charge distribution near the 
charged wall in a system exhibiting critical fluctuations of concentration, remains an open 
question. 



B. Mean-field approximation and continuous Landau-type model 

In the mean-field (MF) approximation the grand potential is assumed to correspond to 
the minimum of the functional: 

fi A/F [ S ,p c ,0,^] = n SR [s,p c ] + (e e [s,p c ,<M - k B TW[s,p c ,cf>}) (17) 

where the electrostatic energy density e e [s, p c , <f>, if)] is given by (jUJ), and k B W[s, p c , 0] is the 
entropy density. In this work we focus on the semi-infinite system, and assume the lattice-gas 
form (ideal mixing entropy) 

w ^ = 1 - p c (z) + S(z) ^ / 1 -p c (z) +S(z) \ + 1 - p c {z) - S(z) k / 1 ~ Pc(z) ~ S(z ) 

2 \ 2 / 2 \ 2 

| p c (z) + (f>(z) ^ ^ p c (z) + (f>(z) j + p c {z) - <j)(z) ^ ^ f) c (z) - (f>( Z ) j ^ 

where we take into account the dependence of the fields on the distance z from the planar 
wall. 

We are interested in thermodynamic conditions corresponding to stability of the uniform 
fluid close to the demixing transition. The equilibrium values of s and p c in the bulk are the 
uniform solutions of the Euler-Lagrange (EL) equations obtained from the minimization of 
the bulk part of the functional (fT7|) ; they are denoted by s and p c . From now on we focus 
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on the deviations from the equilibrium fields 



= s(z) - s (19) 
rj(z) = Pc {z) - p c . (20) 

On the MF level it is possible to take into account the dependence of the dielectric 
constant on the composition and we assume that e given by (jlip depends on the distance z 
from the planar wall according to 

e(z)=e + Se(z), (21) 

where 

<Je(x) = e*0(x) - e„77(x), (22) 

and we introduced the notation 

«. = ^ (23) 

For very small fields i? and 77, the second term in Eq. (|2ip is negligible compared to the first 
one. We shall keep this term in general formulas, valid for arbitrary deviations from bulk 
equilibrium densities. 

The continuum-space Landau-Ginzburg (LG) functional for the fields $(z) and rj(z) is 
defined as 

AC[0, 77, 0, V] = ^ Mi? [s + t?, p c + 77, 0, V>] - tt MF {s, p c , 0, 0], (25) 
where A is the area of the confining surface in the semi-infinite geometry. £[$, rj, <f), ip] can 



be derived from 
in detail in Ref. 



die lattice model in a standard way. Here we follow the method desribed 
. In short, in the first step EL equations on the lattice are derived, 



which contain contributions of the form V 2 / = f(z + 1) + f(z — 1) — 2/(z) for / = s,p c . 
For the boundary layer the effect of the missing neighbors is taken into account, and in 
this way surface EL equations are obtained. At the boundary layer both the bulk and 
the surface EL equations must be satisfied, and the difference of the two equations gives 
the corresponding boundary condition. From the continuous version of the EL equations 
(V 2 / —> d 2 f/dz 2 ) we derive the corresponding functional with the appropriate surface term 
that leads to the boundary conditions. We require that the boundary conditions contain 
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lower-order derivatives than the bulk equations. The functional obtained in this way has 
the form 

C = l ^|^ vT ( 2 ) C ° v W + ^ v,T W Jv/ W + e e + ^Tp| (26) 

+ ^^-M0)w(0), 

where boldface capital letters denote matrices. The transpose of the columnar vectors v and 
v' are v T (z) = ($(z), r](z)) and v' T (z) = (d$(z)/dz, drj(z)/ dz) respectively. The elements of 
the matrix C° = (C,°)i,i=i,2; where indices 1 and 2 correspond to s and p, respectively, are 
given by 



C° ss = k B T 1 P 2 C - QJ SS (27) 



C = k B T[ (1 'J c _ _ 2 QJ PP (28) 



C% = C% = k B T (l _ _' )2 _ _ 2 - 6J ps . (29) 

The second matrix in (126]) is J := (Jy)ij=i,2) again with indices 1 and 2 corresponding to s 
and p, respectively. The electrostatic energy- density e e in the case of a fluctuating dielectric 
constant has the form OH]). The linear surface fields h = (h s , h p ) are 

h s h s JsgS JspPci hp hp JppPc Jps^i (30) 

and P is obtained from the expansion of W about s = s and p c = p c with the linear and 
quadratic parts in the fields $ and rj subtracted. The length unit is the lattice constant 

1/3 

v , comparable with the molecular diameter. The last term describes the electrostatic 
interaction with the uniformly charged surface possessing the surface charge density ea. 



C. Spinodal surface and bulk correlation functions 

In order to find the bulk correlation functions in the Gaussian approximation as well as 
the spinodal surface and the critical line of this model in MF, let us consider the Gaussian 
part of the LG functional in the bulk. In the Fourier representation we have, 

Cg = I Wfl r ^ k) ^ k) (31) 
11 



where v*(k) = v T (— k), and 

C(fc)=C° + it 2 J = J(r 1 C + fc 2 I). (32) 

The correlation functions for the fluctuations of the solvent composition d and the density 
of ions rj in the Fourier representation are given by the matrix G = C _1 . Each component 
of this matrix is inversely proportional to det(J~ 1 C° + k 2 T) = (Af + k 2 )(X 2 + k 2 ), where X 2 
and X 2 , are the two eigenvalues of the matrix 

M = J^C (33) 

and are given by 

TrM =f \l (TrM) 2 — 4 det M 
AL = y{ 2 • (34) 

The correlations are thus given by a linear combination of the two terms, (X 2 + A; 2 ) -1 and 
(A 2 , + k 2 )" 1 . The asymptotic decay of correlations in real space is described by the inverse 
correlation length = min(Ai, A 2 ). 

The uniform phase is stable in the range of parameters T, s and p c such that det C(k) > 
for all k < k max , where k max is an upper cutoff on wave numbers |k| < k rnax = n/a. a is 
usually identified with some appropriate microscopic length, e.g., the lattice spacing or the 
molecular diameter. Because det C(k) = J det (M + k 2 I) = J(A 2 + k 2 )(Xl + k 2 ) where 

J = det J = J ss Jpp - J 2 p , (35) 

for the instability analysis we have to distinguish cases of positive and negative J. 

If J > 0, from the condition detC(O) > (one-phase stability) we obtain after some 
algebra that TrM > and X 2 ,X?, > 0. In that case the instability of the uniform phase 
occurs at k — when 

det C(0) = J det M = JX\X 2 2 = or C ss = (36) 

The condition C ss = determines the instability for the case of the vanishing field t] (the 
same concentration of ions in the phase-separated solvent). From (1341) it follows that X 2 — > + 
so that the asymptotic inverse decay length is = Ai. 

For J < 0, on the other hand, the stability condition of the uniform phase, det C(k) > 0, 
requires det (M + k 2 l) = (Xf + /^(A 2 , + k 2 ) < for all k < k max . In particular, the 
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uniform phase is stable against macroscopic phase separation as long as A^A^ < 0- The 
model can describe phase separation into two uniform phases only when the boundary 
of stability of the uniform phase is associated with k — 0, i.e., when upon decreasing 
temperature (X 2 + k 2 )(\% + k 2 ) changes sign first at k — 0. In our case this means that 
Af A| > (X 2 + fc^ oa ,)(A| + the necessary condition for the above is Af + Ag < 0. This 

condition can only be satisfied when TrM < 0, since in this case A^ < and A2 — > + ; 
hence, 1/C.t = A2. Our model permits an instability at < k < k max , which will lead to the 
modulated phase, for a set of parameters satisfying J < and TrM > 0, since in this case 
A^ — > 0~ and X 2 . > 0. We will not consider such parameters in the current paper. 

For a given composition, the spinodal surface T s (s,p c ) is determined by the higher tem- 
perature satisfying the condition ( |36|) . The explicit relation between the temperature and 
the composition of the mixture that satisfies the first condition in Eq. f[3T)j) . reads 

{k B T s f + k B T s [\2J sp sp c - 6J SS (1 - p c - s 2 ) - 6J pp p c (l - p c )] (37) 
+ 36Jp c ((l - p c ) 2 - s 2 ) = 

For p c = we have simply k B T c = QJ SS and s c = in this model. For a very small amount 
of ions the shift of the critical temperature can be estimated from the derivative of T s with 
respect to p c at s = p c = 0. The derivative at that point is 

dT s Q(J 2 P - J 2 S , ^ 



S'fs ) 



(38) 



dpc J ss 

while the derivative at s = p c = of k B T satisfying C ss (0) = is — 6J SS . Hence, Eq. 
( |37|) describes the actual instability at least for small amount of the solute. The critical 
temperature increases upon addition of ions when J sp — J ss > 0. Experimentally addition of 



some salt leads to the increase of the temperature of the upper critical point [10|, Thus, 
if we require that this model describes the experimental situation on the MF level, we should 
assume that J sp — J ss > 0. This condition, together with the earlier assumption J ss > J pp , 
leads to J < J 2 S — J 2 < 0. On the other hand, for J > the critical point decreases upon 
addition of salt, which is not consistent with the experimental observation. However, the 
present study is limited to MF approximation, and the effect of fluctuations on the critical 
point should be taken into account to draw definite conclusions. We shall consider both 
cases, J > and J < 0, bearing in mind that in the latter case the parameters must satisfy 
the condition TrM < 0. 
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D. Euler-Lagrange equations in semi-infinite system 



Our goal is to determine the shape of the ion number and charge density profiles near the 
charged wall, therefore we need to consider the Euler-Lagrange (EL) equations corresponding 
to the minimum of the functional fl26l) . 

From 5C/5(j)(z) = we obtain 



<p(z) = -p c (z)g(z), 



(39) 



where 



g(z) = tanh(/3eV>0z)/2). 



(40) 



Note that —1 < g(z) < 1. From the above we can obtain ip as a Taylor series in <fi and rj, of 
odd orders in 0. The charge density can be eliminated by using (13"9"1) . and after some algebra 
we obtain three EL equations for rj and g: 



J(ti" + 6tf) 



(J PP + </ ps ) hi(l -p c + s) + {J ps - J pp ) ln(l -p c -s) (41) 



-2 J ps ln(p c ) - J ps ln(l - 5- 2 ) - J pp Ap + J sp Ap c 



J(p" + 6p) 



knT 



(J ps + J ss ) ln(l - p c + s) + ( J ss - J ps ) ln(l -pc-s) (42) 



-2J SS ln(p c ) - J ss ln(l - g 2 ) + J sp Ap - J ss Ap c 



g [e + - e v ri) + g ( e + e ^ ~ e ^ 



2vre 



k B T 



l-g 2 )g 2 Pc = 0, (43) 



where 



Ap = A/i + 6 J ss s + 6 J ps p c 



In 



1 - Pc + S 



1 - p c - S 



(44) 



and 



Ap c = Ap c + 6 J pp p c + 6 J ps s 



k B T 



In 



2 V{l-p c f-p 



(45) 
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In the above equations we used the explicit lattice-gas expressions for the entropy (ITS]) . 
Equations for A/2 and A/2 C follow from the bulk EL equations. 

The boundary conditions of the EL equations for the functional (126]) are 

e + e#&(0) -6^(0) 

— Vip{z)\ z=0 = -ea (46) 

and 

d'-$ = H#, v'~V = H v , (47) 

where 

hpJp S h s Jpp _ hpJp S h s Jpp h s Jp S hpjgg _ h s Jp S hpJ ss 
H* = = s + , H, = = p c + . 

(48) 

J is defined in ( 155]) . 

When e — e, the third EL simplifies, 

g"(z) - 2gg(z)(e/k B T) + {2ne/e){e/k B Tfp c {z)g{zf{l - g{zf) = 0. (49) 

Alternatively, ip could be eliminated, and three EL equations for the fields rj and <fi would 
be obtained. 



III. APPROXIMATE SOLUTIONS OF THE EL EQUATIONS 

The nonlinear EL equations (I41|) - ()47|) can be solved numerically. Since in the critical 
region and for weak surface fields and small a the magnitudes of the fields 0, d and rj are 
small, we can obtain approximate analytical solutions within the perturbation method. 
Analytical solutions can give more general insight, but in the future studies should be 
supplemented with numerical solutions for particular choices of model parameters. 

We postulate that the solutions of the EL equations can be written in the form 

N 
n=l 

where / = 77, (ft, ip, = 0(u n ), and v is a small parameter. 
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A. Solution of the linearized EL equations 



Let us first consider the linearized equations (j4ip - (j43p . For the linearized equations we 
shall simplify the notation, $W = rf 1 ' = 77, = 0. For the electrostatic part we obtain 
the linearized equations of the well known form (the spatial dependence of permittivity leads 
to nonlinear contributions to the EL equations) 



ld2 ' Hz Kem = (51) 



4:Ti dz 2 



which together give 



e^(z) + k -^l = (52) 

Pc 



<j>{z) =K 2 <f>(z), (53) 



where 

Note that n is the inverse Debye length in units of the molecular size. Solution of Eq. (f53"|) 
with the boundary condition Eq. (1461) has the well known form 

(j)(z) = -Kae~ KZ . (55) 

The equations for (rj, $) are formally the same as in the Landau theory for a mixture near 
the demixing critical point, and can be written in the form 

v" = Mv (56) 

where the matrix M = (M^) with the indices 1 and 2 corresponding to s and p respectively, 
is defined in (1331) . and the vector v is defined below Eq. (1261) . In the semi- infinite systems in 
the one-phase state $(z), r](z) — > for z — > 00, and from the Ansatz: 

$ = t\ exp(— \iz) + t 2 exp(— \ 2 z) (57) 
rj = rii exp(— Aiz) + n 2 exp(— X 2 z) (58) 
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we obtain 



(59) 



m M 21 

For T — > T c we have 1/^ — > where = Ai or 1/^ = A 2 for J > or J < respectively 
(see sec JII CI) , and the equations fl5T|) take the form 

0(z)~i ie xp(-z/&) (60) 
?j(z)~niexp(-z/^). (61) 



From the boundary conditions and Eq. fl59|) we obtain in this case 

ni ^t c ^ + ... (62) 

M 22 

h -t-.t c -nijj- + (63) 

M 2 \ 

Note that the linearized equations for <fi and (77, d) are decoupled. Thus, at the linear 
order there is no effect of the concentration profile on the charge distribution. Similar result 



was obtained recently in Ref. 19j. Note, however that we consider a simplified model of SR 
interactions, Eq.flS]), where direct SR couplings between the charge and the concentration 
are disregarded. 



B. Leading-order correction to the linearized EL equations 

Beyond the linear order the coupling between the EL equations leads to modifications of 
the charge and concentration profiles. Note that the solution of the nonlinear equations can 
be written in the form (150|) when all amplitudes of the solutions of the linearized equations 
are small and of the same order, na = 0(u), = 0{y) and U = 0(u). This means a large 
Debye screening length and a small surface-charge density, as well as weak surface fields 
H V ,H S (see Eqs. (|6"2l (I6"31 ). We shall limit ourselves to such conditions and calculate the 
first correction $( 2 \ r/ 2 \ to the solutions of the linearized equations. 

For the electrostatic part we obtain 

(64) 



(2) 



k B T 

epc 



07/ 
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and 



(0( 2 ))" + K 2 ^+^(z) = 



(65) 



where 



H 



{ej&'iz) - e v r)' (z))i/)' (z) + {e#&{z) - e v r](z))ip" (z) 



{(f>{z)r,(z))" (66) 



Note that at the second order in the perturbation expansion we need to calculate the function 
Mcf, to the quadratic order in v. The first term in Eg. (166]) is of the order 0(z/ 2 k 2 ) (recall that 
we limit ourselves to boundary conditions such that the fields rj, <fi obtained from linearized 
equations are all of the order of 0{v)). The perturbation expansion (1501) is justified when 
kg = 0(u) (see (1551) ). i.e., for small surface charge and/or weak screening. When the 
additional condition, namely k 2 <C 1 is satisfied (weak screening), the first term in (16"6"|) 
is negligible compared to the second one. In this case of weak screening the position- 
dependence of the dielectric constant does not play a dominant role for the leading-order 
correction to the solution of the linearized EL equations, and we may assume e = e. From the 
forms of the fields r], <j) obtained from linearized equations we obtain the explicit expression 
for Ncj,{z), and Eq. (l6B1) simplifies to 



) + 



K 



2 ^( 2 ) 



Kae 



(1/6 + n) 2 n ie 



The solution is 



0. 



(67) 



S(2), 



z) = A e~ KZ + A£ b e 



(68) 



where 



A 



(69) 



(C 1 + 2*0 

and .4.0 can be determined from the electroneutrality condition, 

dz(<f) ( - 1 \z)+<J) ( - 2 \z))+a = 0. (70) 

The above condition is satisfied by the solution ( )55|) of the linearized equation, hence for 
the correction term we have the condition J °° dz(j)( 2 \z) = that yields 

Mb 



Aq = —K- 



(71) 
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In the critical region £ b 1 — > 0, hence A£b oo. However, <p^ 2 \z) can be written in the 
following form 

^> M = Ae- + + f e-^f^l (72) 

where 

1 + 

Thus, <p( 2 \z) is regular for ^ b ~ 1 — >■ 0, but the second term in (1721) is non-monotonic as a 
function of z. The ratio between the correction and the leading-order term has the form 

0(D (z) ~ 1 + 

where 



(74) 



l-(l + 6K)(l-e-*/ & ) 



(75) 



is independent of the amplitude n 1; i.e. (f>^ (z) / (j)^ (z) oc rii. The coupling between the 
concentration and the charge fluctuations in the functional (12"6T) leads to the non-monotonic 
correction to the linear solution of the EL equations, with the extremum at 

;W = 2&ln^±^\ (76) 

From Eq. (]74p it follows that the correction term changes its sign at z = z extT /2, i.e., the 
critical adsorption for z < Zq leads to the charge density larger (smaller) than that predicted 
by the Debye-Huckel theory for ri\ > [n\ < 0), whereas for z > zq the charge density is 
smaller (larger), respectively. The effect is proportional to the amplitude n\ of the deviation 
of the density of ions from the bulk value, p c {z) — p c ~ ^i exp(— z/^). Figs. [T] and [2] show 
examples for the specific choices of the model parameters close to the critical point of the 
binary solvent, i.e., for 3> 1. The rescaled ratio R between the leading-order correction 
and the Debye-Hiickel approximation is shown in Fig. [3] for four apirs of the parameters 
(«,&)■ 

Let us consider the non-monotonic correction term for k <C ^ b ~ 1 , i.e., away from the 
critical point. From (]74|) it follows that it has a simple exponential form with the amplitude 
equal to the amplitude ny. 

(f)^(z)/(j) w (z)^n l e- z/ ^ «<C 1 - ( 77 ) 
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FIG. 1: (Color online) Charge density profiles in the semi-infinite system. The dashed line is the 
solution of the linearized equation (|55p . (j)^(z), the dash-dotted line is the correction term (|74|) . 
<fi( 2 \z), and the solid line is the sum of the two functions, <f>^ l \z) + (f)( 2 \z), for k = 0.01, a = 0.5, 
£b = 100, and n\ = 0.5. The charge density is in e/vo units, where e is the elementary charge, and 

1/3 

z is in v Q units, where vq is the volume per molecule in the close-packed system and we assume 
that all molecules are of similar size. 

For T — > T c we have — > 0, and using <*twe obtain the approximation 



^\z)/^\z) ~ | 



1 - « 



(1 - e - */ 6 



« » c 1 



For 2 ^ the above takes the simple form 

^ 2 \z)/(f)^(z) ~ 
The correction term changes sign for zq ~ l//c. 



ni 
2 



1 — 



(78) 



(79) 



Equations for the first corrections to the nonelectrostatic part read 

(v( 2 ')" = Mv< 2 »+D, (80) 
where v^ 2 ^ = ($( 2 \ rj^), and the components of the vector D T = (D s ,D p ) are 

u s - J (81) 

n D b ^ 2 + D'JV + D b m V 2 + D\^ 2 

U p - - (82J 
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FIG. 2: (Color online) Charge density profiles in the semi-infinite system. The dashed line is the 
solution of the linearized equation (|55p . <p^(z), the dash-dotted line is the correction term (|74|) . 
(j)^(z), and the solid line is the sum of the two functions 4>( 2 \z), for k = 0.01, a = 0.5, = 100, 
and n\ = —0.5. The charge density is in e/vQ units where e is the elementary charge, and z is in 
v^ 3 units, where vo is the volume per molecule in the close-packed system and we assume that all 
molecules are of similar size. 
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FIG. 3: (Color online) The rescaled ratio R (Eq. (|75j) ) between the leading-order correction Q74|) 
and the Debye-Hiickel form (|55p of the charge density profile for several choices of pairs («, as 

1/3 

indicated in the plot, z is in v Q units, where vq is the volume per molecule in the close-packed 
system and we assume that all molecules are of similar size. 
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As follows from Eq. (]80[) the correction terms to the composition profile $(z) of the 
mixture are of the order of e~ 2z ^ b and e~ 2KZ . Note that the effect of the charges on the 
profiles of the composition of the mixture is negligible in the critical region, where k ^> 

IV. SUMMARY AND OUTLOOK 

Starting from the lattice-gas model of a four-component mixture we have derived the 
continuum Landau- Ginzburg model for binary mixture solvents in the presence of ions near 
the critital point of the demixing transition. The model encompasses the composition of 
the binary solvent field the density of ions field rj and the charge density field </>. It 
takes into account electrostatic interactions and the preferential solvation. The coupling 
constants appearing in this extended Landau-Ginzburg theory are given explicitly in terms 
of thermodynamic quantities, the temperature, the mean composition of solvent, the mean 
density of ions, and the interaction parameters Jy characterizing the lattice-gas model of a 
mixture. We have assumed that ions are of similar chemical nature. 

The main difference between our functional and the functional studied in Ref.j^] is the 
presence of the term oc rf and terms oc (V77) 2 and oc Vr/V^ in fl26|) . which result from 
short-range interaction potentials. These terms lead to the mixing of the fields $ and rj in 
the critical order parameter. In the semi-infinite system these terms are important for the 
form of the profile of the field 77, and through the coupling of 77 with the charge density </> in 
the entropy term, they influence the charge profile <j)(z). In our approach, direct couplings 
between the charge and the concentration are disregarded, but such coupling would be 
present in the case of ions of different chemical nature (compare Eqs.([5]) and ([8])). 

Mean-field theory for our Landau-Ginzburg model yields the shift of the critical point 
of the demixing transition with respect to the case of binary solvents without ions. The 
direction of the shift depends on the relative strength of the ions-solvent and solvent-solvent 
interaction parameters J sp and J ss , and is positive for J sp > J ss . The linearized EL equations 
in the presence of a charged wall do not lead to the effect of the concentration profiles on the 
charge distribution. We treat nonlinear effects using a perturbation expansion, which gives 
the simple expression for the leading correction to the solution of the linearized EL Eqs. (P74]) 
and fT75|) . The ratio between the next-to- leading and leading terms is proportional to the 
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amplitude n\ of the decay of the ion density, 77 = n x exp(— z/£&), and otherwise depends 
only on the dominant lengths in the system, the bulk correlation length and the Debye 
screening length 1/k. We find that due to the critical adsorption of that component of the 
solvent in which the ions are preferentially dissolved, the amount of counterions in the layer 
near the charged surface of the thickness of the Debye screening length is increased with 
respect to the one away from the critical point. Critical adsorption of the better solvent 
enhances screening whereas the critical adsorption of the poorer solvent leads to the opposite 
effect. This approach is quantitatively valid provided the Debye screening length is large, 
the surface- charge density is small, and the surface fields are weak. We expect the same 
trend beyond the approximate perturbation expansion solution. This will be studied further 
by the numerical treatment of the full EL equations. 

We should note that our theory shows that the leading-order correction to the Debye- 
Hiickel form of the charge profile is determined by the average dielectric constant, i.e., the 
spatial variation of the dielectric constant can be neglected when k <C 1 (in units of inverse 
molecular size). 

Finally, the change of the charge distribution near the charged surface occuring upon ap- 
proaching the critical point of the binary solvent, indicates that in the confining geometries, 
i.e., for systems between two surfaces, the effective interactions between the confining sur- 
faces will be altered. Our predictions will be tested in the future work by explicit calculations 
and by comparison with the experimental data of Ref. \t\. 
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